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Abstract 

The origin of microinertia of micromorphic theories is investigated from the 
point of view of non-equilibrium thermodynamics. In the framework of dual in¬ 
ternal variables microinertia stems from a thermodynamic equation of state related 
to the internal variable with the properties of mechanical momentum. 


1 Introduction 

The prediction of a thermomechanical behavior of solids with microstructure cannot 
be completely correct unless the influence of the microstructure is taken into account. 
In spite of the possible small scale characterizing the microstructure, its influence may 
be not necessary small if corresponding gradients are comparable with those at the 
macroscale. It is well understood that microstructural effects can be described only by 
means of a generalization of classical thermomechanical theory. Such a generalization 
is usually achieved at the expense of the extension of the state space by additional 
independent variables like the microrotation in Cosserats media, the microdeformation 
in micromorphic media, or the heat flux in extended thermodynamics. 

Making the distinction between the medium and its microstructural components, 
one implicitly admits that they can respond to dynamic loads differently. This is re¬ 
flected in the notion of microinertia appeared in the description of micromorphic me¬ 
dia Eli. Mindlin [2j] described the deformation of what he called "cell" explicitly, 
while Eringen and Suhubi |U] introduced an averaging procedure over the microele¬ 
ment. Both approaches results in similar balances of micromomentum in the linear 
case. As noted by Eringen and Suhubi [l|], the corresponding microinertia tensor "re¬ 
sembles" the physical inertia. This means that the coefficient of the second time deriva¬ 
tive of the microdeformation is consistent with the constitutive relations. 
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Mindlin as well as Eringen and Suhubi considered their microdeformation as an ad¬ 
ditional degree of freedom. Therefore, it should be governed (controlled) by boundary 
conditions pointed out in the general form by both approaches. However, it is difficult 
to imagine boundary conditions that excite the microstructure only, without affecting 
the medium itself. This difficulty brings the internal degrees of freedom together with 
internal variables of state, which even cannot be controlled by boundary conditions. 

The classical treatment of internal variables suggests that their evolution equations 
may contain only the first-order time derivatives, and, therefore, have no relations to 
inertia [3r[6]. This is definitely true for order parameters considered as internal vari¬ 
ables, and for pure dissipative effects like heat conduction. However, the general¬ 
ization of the internal variables theory by the dual internal variables concept allows 
hyperbolic evolution equations as well IZQ. The question is how the coefficient of the 
second order time derivative of the microdeformation corresponds to microinertia. 

In order to answer this question, we will show how the formalism of dual internal 
variables works. Starting with simplest possible case, we demonstrate the difference 
between the standard single internal variable theory and the dual internal variables 
concept. Then we compare the generalization of linear elasticity by dual internal vari¬ 
ables with the Mindlin microelasticity model [3]. For simplicity, all the considerations 
are made in one-dimensional case. 


2 Single internal variable: one-dimensional example 

While for mechanical motion of a material point it is sufficient to prescribe its coordi¬ 
nates and velocities, the theory of internal variables is essentially thermodynamical. 
This means that we deal with a thermodynamical system and suppose that all thermo¬ 
dynamic quantities like temperature, energy, entropy, etc. are defined. 

2.1 Entropic representation 

Here we assume that the entropy density S is specified as a function of the internal 
energy density E and an internal variable cp and its space derivative (p x - 

S = S(E,cp,cp x ). (1) 

Constitutive assumption ([T|) allows us to introduce the corresponding conjugate quan¬ 
tities as follows: 

1 _ as r _ as rj _ as 

9'~dE' 6 dq>' 9 dy~ x ' ( } 

where 6 is the temperature. With the usual thermodynamic notation with differentials 
the consequent Gibbs relation can be written as 

dE = 9dS — rdcp — rjd<p x . (3) 

Therefore, Eq. (QJ is the equation of state of a "gradient" (or weakly non-local) internal 
variable theory. 
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The conservation of the internal energy can be represented as 


Et + Qx — 0 , 


(4) 


where Q is the heat flux, index t denotes time derivative. For the calculation of the en¬ 
tropy production a generalized entropy current density (i.e. entropy flux) is introduced 
in the form 


/ = 


Q - yq> 

e ' 


(5) 


Here the deviation of the entropy current density from its classical form J c i ass = Q/6 
is represented by a particular additional term. This form is characteristic for second- 
order weakly nonlocal theories and can be justified either by a rigorous exploitation of 
the entropy inequality [8,@], or by the classical separation of divergences [10, 6]. In the 
next section we will demonstrate the second method. 

With the help of modified entropy flux © the entropy production can be calculated 
as follows: 


St + Jx 


1 _ T 11 

-E t + -cp t + -Q<Ptx + 


Q-yq> 

e 




y_ 

16 


(Q-y?t) ( ^ +j(t-t] x )> o. 


( 6 ) 


A constitutive relation for heat flux and an evolution equation for internal variable 
are needed to be specified for any considered process. The given form of the entropy 
production supposes two possible representations of the thermodynamic fluxes and 
forces as one can see inspecting the last two lines of Eq. ©. In the first form the thermal 
part is the classical, while in the second form the thermal part is modified. However, 
the linear constitutive equations are identical, as it is shown in the Appendix. 


2.2 Energetic representation 

Free energy is the Legendre transformation of the internal energy by the entropy W = 
E — SO. It is convenient for practical use, because temperature is the introduced vari¬ 
able of state instead of entropy or internal energy: 


W = W(6, f/fx ). 


Partial derivatives of W define 

c ._ 

^ ~W' T ~~ 9^' ?/ 


(7) 

( 8 ) 




which results in the Pfaffian form 
dW = — Sd9 — rdcp — i]dq) x . 


(9) 











The energy conservation equation is 

(S6)t + Q X = h int , h int := -Wt, 


( 10 ) 


where the right-hand side of Eq. (TlOlh is formally an internal heat source 

The energy conservation equation should be accompanied by the second law of 
thermodynamics, here written as 


St + {Q/e + K)x> 0, 


( 11 ) 


where in the entropy flux is additively separated onto the classical part and the extra 
entropy flux K. Multiplying the second law ([Til by 8 

8S t + 8(Q/8)+K) x >0, (12) 

and taking into account Eq. (fTOb . we obtain 

- (W t + S8 t ) + (6K) X - ( Q/8 ) + K)8 X > 0. (13) 

The last equation can be represented in the form 

S8 t + ( Q/8 ) + K)8 X < h int + ( 8K) X . (14) 

The internal heat source h int is calculated following the constitutive assumption (0 


h int = -W t = 


dW n dW dW 

-8t - —<pt 


cp xt = S6 t + r ft + V?xt = h th + h 


Jntr 


d6 d(p dcp x 
Accounting for Eq. (fT5l) , dissipation inequality (IT4l) can be rewritten as 

O = rcpt + tj(p x t — (Q/ 8 + K)8 X + ( 8K) X > 0. 


(15) 


(16) 


Here O is the entropy production multiplied by the temperature. To rearrange the 
dissipation inequality, we add and subtract the same term rj x (pt 


O = X(p t + ij(p x - rj x cp t + r/ x cpt - ( Q/8 + K)8 X + (8K) X > 0, 
that leads to 

= (t - tj x )(pt - ( Q/8 + K)8 X + ( tjcpt + 8K) X > 0. 


(17) 


(18) 


Therefore we can define the "extra" entropy flux by the elimination of the divergence 
term in Eq. (IT8l) 


K= -8 1 ycpt- 


(19) 


This method, the separation of divergences, has its roots in the classical irreversible 
thermodynamics [10] and formulated explicitly in case of internal variables by Maugin 
[H]. Then the dissipation inequality reduces to 


= (t - y x )(pt - (Q + ycpt)(log(6)) x > 0. 


( 20 ) 


It is remarkable that in the isothermal case (8 X = 0) the dissipation is determined by 
the internal variable only. It is easy to check that inequalities (l20l) and ((6]) are identical. 
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2.3 Evolution equation for a single internal variable 

In the isothermal case the dissipation inequality (1201) is 

®=(T-t] x )<p t >0, (21) 

The linear solution is 


cpt = k (r - rj x ), k> 0, 

since dissipation inequality (1211) is satisfied automatically in this case 
<E> = kcp^ >0, if k > 0. 


( 22 ) 


(23) 


It is easy to see that the dissipation is the product of the thermodynamic flux cp t and the 
thermodynamic force (r — ;/ T ). The proportionality between the thermodynamic flux 
and the conjugated force is the standard choice to satisfy the dissipation inequality 
To see how the obtained evolution equation looks like, we specialize free energy 
dependence CD) in the isothermal case to a quadratic one 

W = 1 Bcp 2 + ±C<pl (24) 

where B and C are material parameters. It follows from equations of state © that 


aw 


r : = 


= -Bcp, rj : = 


aw 


— C(pXr 


dcp T ' ' • dcp , 

and evolution equation (122)) is an equation of reaction-diffusion type 
cpt = k(Ccp xx - Bcp), k> 0. 


(25) 


(26) 


The given standard formalism of the single internal variable of state is sufficient for 
many cases [12-19]. In the case of a more general free energy, Eq. (1261) is the Ginzburg- 
Landau equation [20(]. 


2.4 Dual internal variables in one dimension 


Now we would like to extend the same technique to the case of two internal variables. 
However, none of them is presupposed to be definitely an internal variable of state or 
an internal degree of freedom. Let us consider them as internal variables "in general". 
Here we treat the isothermal case, therefore we extend the free energy focused train of 
thought of the previous section. The corresponding entropic treatment is straightfor¬ 
ward and equivalent. 

Let us suppose that the free energy density depends on the internal variables cp, ip 
and their space derivatives 

W = W(8, cp, cp x , ip, ip x ). (27) 


The equations of state in the case of two internal variables read 


aw 

He' 


aw __ aw aw 

lip' >] ~dy~ x ' ^ ~^p' ^ 


aw 
a< px 


r : = 


(28) 



















The Gibbs relation follows accordingly 

dW = — Sd6 — rdcp — i]d(p x — £dtp — £ dip x . 


( 29 ) 


Then we consider a non-zero extra entropy flux following the case of a single internal 
variable and set 


K = -e- 1 tj( Pt -e- 1 tf t . 


(30) 


It can be checked that the intrinsic heat source is determined in the considered case as 
follows 


h intr := (t -Tj x )(p t + ^-^ x )xpf 


(31) 


The latter means that the dissipation inequality in the isothermal case reduces to 

= (t - rj x )(pt + (£- £ x )ipt > 0. (32) 

The solution of the dissipation inequality can be represented as [@] 


T)= iTl-’h, where L = T 1 ^ 


Ipt 




L 2 1 L 2 2 


(33) 


Nonnegativity of the entropy production (132J) results in the positive semidefiniteness 
of the symmetric part of the conductivity matrix L, which requires 


L n >0, L 22 > 0, L n L 22 - (Ll2 + L2l)2 >0. 

To be more specific, we keep a quadratic free energy density 
W = iBy 2 + lapl + 1 Dip 2 + jFipl 
Calculating the quantities defined in Eq. 


r := 


£ ; = 


aw 

dcp 

aw 


= -B<p, 7 / : = 


= -Dip, l := 


aw 

dcp x 

aw 


— Ccpx, 

= —F ipx, 


dip dip x 

we can represent system of Eqs. (l33l) in the form 

cpt = Ln( — Bcp + Ccpxx) + ki 2 ( — Dip + Fip X x), 
ipt = h 2 i (—Bcp + Ccpxx) + h 22 (—Dip + Fip X x). 


(34) 

(35) 

(36) 

(37) 


(38) 

(39) 


Now we will derive a single evolution equation for the primary internal variable. To 
this end a suitable rearrangement of the previous equations is performed: 


l i + Ln(B-cFL)) V + L 12 [D-F 

3 2 


MB-Cj^k+I^ + MD-F 


ax 2 
a 2 , 


dt 


dx 2 


ip = L n <p + Li 2 ip = 0, (40) 

ip = L 2 i<p + L 22 i/7 = 0. (41) 


A 











Here the traditional notation of time and space derivatives is used for the clear dis¬ 
tinction of the differential operators L\\, L\ 2 , L 21 and L 22 . Then the easiest way to 
eliminate xp is to multiply the first equation, (1401) , by L 22 and the second one, (l4l|) , by 
L \2 and subtract them: 


L 22 L 11 CP - L 12 L 21 CP — (pu — detL(CD + BF)cp xx + 

+ (BLn + DL 22 )cpt — ( CLu + BL 22 )(ptxx + det L ( BDcp + CFcp xxxx ) = 0, 


(42) 


where detL = Ln L 22 — I 42 T 21 is the determinant of L. The free energy density W 
is non-negative by default, which results in non-negativity of material parameters 
B, C, D, and F. This means that Eq. (l42l) is the hyperbolic wave equation with dis¬ 
persion and dissipation. 

Thus, extending the state space of our thermodynamic system by an additional in¬ 
ternal variable and keeping the quadratic form for the free energy density, we arrive at 
the hyperbolic evolution equation for the primary internal variable. The correspond¬ 
ing evolution equation for the secondary internal variable can be derived similarly. 
The natural question is the following: can we associate the abstract internal variables 
with the description of any physical process? 

2.5 Example I: Linear elasticity 

To answer the question posed above, we simplify the situation as much as possible. We 
note first that dissipation will be absent if the diagonal components of the conductivity 
matrix are zero (Ln = L 22 = 0). Due to inequalities (l34l) this requires antisymmet¬ 
ric conductivity matrix, that is, Casimir type reciprocal relations L \2 = — L 21 . Next, 
dispersion will be eliminated if we choose the material parameters B = 0 and F = 0. 

It follows immediately that the evolution equations for the dual internal variables 
(l38l) and (l39l) can be rewritten as 


(43) 


cpt = -L 12 Dxp, 


ipt = -Li 2 C(<p x ) x . 


(44) 


Without loss of generality, we can choose the value L 12 = — 1. We should emphasize 
that the two evolution equations express the duality between internal variables: one 
internal variable is driven by another one and vice versa. 

Now we are ready to relate the introduced internal variables to well-known physi¬ 
cal quantities. First, we note that Eqs. (l43l) . (l44l) can be represented in the form 


cp xt = Dip X/ 


(45) 


ipt = C{cp x ) x . 


(46) 


Secondly, denoting cp x as e, and ip as pv we arrive at the system of equations 

£t = D(pv) x , 


(47) 


(pv) t = Ce x , 


(48) 
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which coincides with one-dimensional equations for linear elasticity (with the standard 
notation for derivatives), and p is constant: 


de dv 
dt dx' 


(49) 


dv 3(7 
P dt = dx' 


(50) 


after the following choice of coefficients: D = 1 /p and C = pc 1 , where c is the elastic 
wave velocity and a = pc 2 e. This means that the internal variable cp can be interpreted 
as the displacement and the internal variable ip as the momentum. These variables are 
controlled by boundary conditions and, therefore, they are true mechanical degrees of 
freedom, while appeared as internal degrees of freedom first. No doubts arise about 
inertia in this case. 

We have considered a thermodynamic system (body) under assumption that there 
may be certain internal phenomena which we tried to describe by internal variables. 
We have derived possible evolution equations for these internal variables from the 
first and second laws of thermodynamics. It has been demonstrated that evolution 
equations for internal variables may be hyperbolic for a quadratic free energy density. 
Moreover, in the non-dissipative and non-dispersive case these evolution equations 
can be identified with equations of linear elasticity. 

It follows that both internal variables of state and internal degrees of freedom can be 
considered as particular cases in the framework of unified formalism. This formalism 
can help us to consider generalized continuum models by means of internal variables. 
In our next example, we look at the Mindlin micromorphic elasticity theory. 


3 Example II: One-dimensional microelasticity 

Following I@], in our one-dimensional treatment we introduce the macrostrain e, the 
microdeformation ip, the relative deformation 7 = e — ip, and the microdeformation 
gradient ip x ■ The corresponding quadratic free energy density [ 0 ] is reduced to 

1 b a 

W(e,y,ip x ) = -pc 2 e 2 + ge 7 + -y 2 + pi . (51) 

Equations of state determine the corresponding stresses 

aw 2 aw aw 

r = — = pc A e + gy, cr = - = ge + by, p = — = aip x , (52) 

and equations of motion have the form [ 2 ] (without body forces) 

putt = T X + <T X , (53) 

1 . ? 

-p d iptt = p x +a, (54) 

where u is the displacement, and p'd 2 /3 plays the role of microinertia. 
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Now we compare the Mindlin microelasticity theory with the dual internal vari¬ 
ables approach. To avoid the confusion with the original notation of Mindlin, we use 
other symbols for internal variables. We suppose that the free energy density depends 
on the gradient of the displacement, u x , the internal variables a, /3 and their space 
derivatives 


W = W(u x ,oc,a x ,p,p x ). 

Then the equations of state follow 

, aw , aw 


aw 


a : = 


r : = 


du x ' doc ' ^ dcc x ' v dp ' * 9/3} 

As in [1|, we use a quadratic function for the free energy density 

_ ^ 111 

W = -pc 2 u x + Acaix + -Boc 2 + -Ccc 2 + -D/3 2 . 


£:= 


aw 


£:= 


aw 


(55) 


(56) 


(57) 


Considering the non-dissipative case, we have for the evolution equation for internal 
variables (cf. Eqs. (138)) , (l39l) with Ln = L 22 = 0) 


— £42 (£ — Cx) — —L 12 D/3, 


(58) 


i &t = -Li 2 (t' - rj x ) = -L 12 (Gr xx - Aw x - Ba). (59) 

As a result, we arrive at the hyperbolic evolution equation equation for the primary 
internal variable 

lot-tt = Coc xx - Au x - Boc, (60) 

where 1 = 1/ (L 2 2 D) is the measure of microinertia. 

The balance of linear momentum 

pv t = <t' x , (61) 

results in the form 

pu tt = pc 2 u xx + Aoc x . (62) 

Let us rewrite the governing equations of the Mindlin microelasticity (l53l) and (f54l) 
using the notation u x instead of e and cl instead of the microdeformation ip. We have, 
respectively, 

putt = T X + (T X = pc 2 u xx + (b + 2g)u xx ~(g + b)cc x , (63) 

'Y 

-p’cfcLtt = a<x xx + (g + b)u x - boc, (64) 

It is easy to see that systems of equations (1601) , (l62l) and (1631) , (1641) coincide if the follow¬ 
ing values of parameters are chosen: A = —(b + g), B — b,C = a,I = 1/ 3p'd 2 , b + 2g = 
0. This means that the dual internal variable theory can recover the generalized contin¬ 
uum theory of Mindlin ( 2 ] and the corresponding microinertia measure has the same 
justification as in the Mindlin microelasticity. 


q 







4 Summary and discussion 


In this short note we have analyzed the relation of mechanical and thermodynamical 
approaches of extending the state space of continuum theories. Our analysis focused 
on the origin of inertial terms in the evolution equations. In the first part of the paper 
we have presented the thermodynamic background of single and dual internal variable 
theories both in the entropic and in the energetic representation [21] . 

The recent methodology of dual internal variables unifies these treatments and we 
have shown that both classical linear elasticity and Mindlin type microelasticity are 
well suited in this general background without a direct reference to mechanical princi¬ 
ples or notions. Our treatment was restricted to one spatial dimension and without a 
mention of ispace-time related conceptual questions. The 3D tensorial repr esentation 
of the application of the dual internal variables approach is given in [22, 231- 

What is the origin of inertia in the light of our approach? The two ingredients were 
the thermodynamic equation of state and the antisymmetric part of the constitutively 
constructed dynamics. The more surprising part is the first: in our treatment iner¬ 
tia and microinertia stem from an equation of state. The coefficients C and D in the 
quadratic free energy function of linear elasticity, (Eq. (1351)1, are connected with the 
density, which is the measure of inertia. In case of microinertia the coefficient D is 
enough, the gradient term of the second variable was not necessary. In our dual inter¬ 
nal variable approach one of the variables plays the role of momentum when compared 
to non-dissipative mechanics, on the other hand it is a state variable from a thermody¬ 
namic point of view. 

This kind of close connection of mechanics and thermodynamics is not too surpris¬ 
ing from the broader perspective of relativistic theories. Therein momentum cannot 
be separated from energy in a covariant treatment, therefore an energy dependent en¬ 
tropy depends naturally on the momentum as well. This is the root of the famous 
problem of the proper Lorentz transformation of the temperature. The first formu¬ 
lation of Planck and Einstein considered energy and momentum as variables of the 
entropy [24-2^]. Then Blanusa and Ott seemingly have eliminated this dependence by 
fixing the thermometer to the body [27t|28fl. However, in a proper covariant treatment 
these approaches are only two sides of the same coin, where the entropy depends on 
energy-momentum [29t |30fl. 

In our nonrelativistic treatment the integration of mechanics and thermodynamics 
looks like more involved than usually taken. From the point of view of the presented 
approach it is remarkable, that Hamiltonian dynamics is a particular case characterized 
by zero dissipation. Moreover, the possibility of a pure variational description can be 
a consequence of ideal dynamics without dissipation 01. 

In summary, we have shown that inertial terms are natural in a thermodynamic the¬ 
ory with dual variables and the conditions of their appearance is well understandable 
in terms of mechanical notions. However, a true space-time related frame independent 
treatment is required to a final clarification of this issue. 
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5 Appendix: Representations of entropy production 

In some cases entropy production ©2 appears in the following form 
E = /iX 1 + / 2 (X 2 -AX 1 ) >0, 


(65) 


where ]\ and J 2 are constitutive quantities, and they are to be determined as functions 
of X\ and X 2 , A is a given function of Xi and X 2 . One may introduce the following 
choice of thermodynamic fluxes and forces: 


Fluxes 

Ji 

h 

Forces 

Xi 

X 2 - AXi 


Table 1. First choice of thermodynamic fluxes and forces 
The linear solution of inequality (l65|) results in: 
h = hX 1 + l 12 (X 2 -AX 1 ) r 

h = / 2 iX 1 + / 2 (x 2 -ax 1 ). 


( 66 ) 

(67) 


Here the inequality requires that the symmetric part of the coefficient matrix is positive 
definite, therefore l\, l 2r l \ 2 and l 2 \ coefficients fulfill the following inequalities: 


h >0, l 2 > 0, and l\l 2 — 


{J12 + l. 


2b 


> 0 . 


( 68 ) 


However, the entropy production may be written in a distinct but equivalent form 

x = (Ji - AJ 1 )X 1 + J 2 X 2 > 0. (69) 

This grouping of the terms suggest a different choice of the thermodynamic forces and 
fluxes: 


Fluxes 

h ~ AJ 2 

h 

Forces 

Xi 

X 2 


Table 1. Second choice of thermodynamic fluxes and forces 

The consequent linear solution introduces the constitutive relations as follows: 

Ji ~ AJ 2 = k\X 1 + k\ 2 X 2 , 
h = ^ 21 X 1 +k 2 X 2l 

where the coefficient matrix must fulfill the following inequalities: 
k\ >0, k 2 > 0, and k\k 2 — —- 12 > 0. 


(70) 

(71) 
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( 72 ) 










The material coefficients l\, l 2 , hi and l 2 \ are not independent of k\, ki, k\ 2 and k 2 \ and 
their relations can be expressed easily: 

f h hi\ _ fh +A(k 12 + k 21 ) + A 2 k 2 
V21 h) V ^21 + Ak 2 

The inverted form is also explicit 

k\ ki 2 \ _ fh ~ A(hi +l 2 \) + A 2 l 2 

k 2 \ k 2 J y J 2 i — Al 2 

As one can check, inequalities (l68l) and (l72l) are equivalent and therefore the constitu¬ 
tive equations (l66]) -(l67l) are equivalent to (l70l)-(l7TI). The seemingly different thermody¬ 
namic fluxes and forces lead to equivalent constitutive relations. 


k\ 2 + Ak 2 
k 2 


hi — Al 2 

h 


(73) 


(74) 
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